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Abstract. We present a first-principles density-functional calculation for the Raman spectra of a neutral BEDT-TTF 
molecule. Our results are in excellent agreement with experimental results. We show that a planar structure is not a stable 
state of a neutral BEDT-TTF molecule. We consider three possible conformations and discuss their relation to disorder in 
these systems. 

Crystals based on the BEDT-TTF molecule of the form (BEDT-TTF^X, where X is a monovalent anion 
such as Cu[N(CN) 2 ]Cl, are of great interest both experimentally and theoretically. This is due to the 
variety of ground states found, and the sensitivity of the ground state to subtle changes in pressure and 
chemistry pQ. It has been proposed that a minimal model for such crystals is the Hubbard model on an 
anisotropic triangular lattice for which we ultimately aim to extract the parameters from quantum 
chemistry calculations, using NRLMOL, a Gaussian based, all electron density functional theory based 
electronic structure method developed by Pederson et al. 

We have calculated the Raman spectra of a neutral BEDT-TTF molecule. We have assumed the 
molecule is symmetric about the x, y and z axes, that is, that it is in a planar structure, shown in figure 
2. The results are shown in figure 1. The range of the plot was chosen to aid comparison with the 
experimental data from Eldridge et al. [I] reproduced in figure 1. Outside the range of frequencies shown 
in figure 1 the only notable features are several strong hydrogen stretch peaks around 3000 cm -1 that were 
also reported by Eldridge et al. [4J. There were also several peaks present with imaginary frequencies, 
which show that the planar structure is not a stable state of the BEDT-TTF molecule. This agrees with 
the results of quantum chemistry calculations by Demiralp and Goddard jSj. 



0.06 



30 r 



25 



£20 
CO 
■* 

<15 
>. 

CO 

§10 




1493 — 



BEDT-TTF 
Room Temp. 1509 



990 1551 
1000 
1013 1422 
1125 1408 
1174 
1257 
i 1283 



500 



1000 



1500 




v cm 



200 400 600 800 1000 1200 1400 1600 
Wave Numbers (cm- 1 ) 



Figure 1: Left: Raman spectra for BEDT-TTF at 300 K calculated using the generalised gradient approximation form 
of density functional theory methods. The peaks are broadened by a Gaussian curve with a full width half maximum of 1 
cm -1 to aid comparison with experiment. Right: experimental Raman spectra for BEDT-TTF at room temperature from 
Eldridge et al. @J. 

As our calculations have shown that the planar structure is not stable, we have performed calculations 
on several other conformations of the BEDT-TTF molecule with no symmetry enforced. Two such struc- 
tures of these are shown in figure 2. The total energy of the two alternative conformations are similar to 
each other at about 1 eV below that of the planar structure. 

It has been shown that conformational changes in the terminal ethylene groups of BEDT-TTF can 
lead to intrinsically non-magnetic disorder [Hill]- In both K-(BEDT-TTF) 2 Cu[N(CN) 2 ]Br and /?-(BEDT- 
TTF)2l3 this disorder has been shown to lead to a suppression of the superconducting critical temperature 
|3 IB] • Understanding this phenomena and the details of the conformational change involved may therefore 
be important in identifying the pairing symmetry of the superconducting state of these materials [7J. It has 
been proposed by Ravy et al. [9J that the conformational changes can be modeled as a double harmonic 
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well. Disorder becomes frozen in below T - 80 K in K-(BEDT-TTF) 2 Cu[N(CN) 2 ]Br and k-(BEDT- 
TTF)2Cu[N(CN)2]C1 jSJ. It has also been found that the activation energy for terminal ethylene disordering 
is ~ 3200 K in /c-(BEDT-TTF) 2 Cu[N(CN) 2 ]Br and ~ 2650 K for k-(BEDT-TTF) 2 Cu[N(CN) 2 ]C1 0. The 
calculated energy difference between the W and Z structures (see figure 2) is 750 K. Clearly the potential 
barrier between the relevant conformations then must be greater than the energy difference between the 
two conformations. However, the activation energy will be determined by the highest energy state which 
the molecule has to pass through to move from one configuration to another. In this case of a BEDT-TTF 
molecule moving between the W and Z configurations this may well be close to the the energy difference 
between the W conformation and the planar conformation, 21200 K, which is an order of magnitude larger 
than the measured activation energy. On the other hand, it is possible that there is path between the 
W and Z configurations that does not involve moving through the planar configuration and has a lower 
activation energy. Also it is clear from the difference between the activation energies of the two salts 
that interactions with the local environment can change the activation energy of conformational changes. 
Therefore, although although it seems unlikely that the W and Z conformations are responsible for terminal 
ethylene disorder in BEDT-TTF salts, this cannot be strictly ruled out. 
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Figure 2: The conformations of the BEDT-TTF molecule considered in this paper and the relative energies and HOMO- 
LUMO gaps. 

In conclusion, we have calculated the Raman spectra of a neutral BEDT-TTF molecule using first- 
principles density functional theory methods. Our calculated spectra agrees well with experimental results. 
The presence of peaks with imaginary frequency show that the structure enforced on the molecule for 
the calculation is not the stable state. We have shown that the energy difference between the W and Z 
conformations is not inconsistent with the reported activation energy of K-(BEDT-TTF)Cu[N(CN) 2 ]Br and 
fi>(BEDT-TTF)Cu[N(CN) 2 ]Cl, but a crude estimate of the activation energy between these conformations 
is an order of magnitude larger than the experimentally observed activation energy. 
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